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The Kondo zero bias anomaly of Co adatoms probed by scanning tunneling microscopy is known 
to depend on the height of the tip above the surface, and this dependence is different on different low 
index Cu surfaces. On the (100) surface, the Kondo temperature first decreases then increases as 
the tip approaches the adatom, while on the (111) surface it is virtually unaffected. These trends are 
captured by combined density functional theory and numerical renormalization group (DFT+NRC) 
calculations. The adatoms are found to be described by an S' = 1 Anderson model on both surfaces, 
and ab initio calculations help identify the symmetry of the active d orbitals. We correctly reproduce 
the Fano lineshape of the zero bias anomaly for Co/Cu(100) in the tunneling regime but not in the 
contact regime, where it is probably dependent on the details of the tip and contact geometry. The 
lineshape for Co/Cu(lll) is presumably affected by the presence of surface states, which are not 
included in onr method. We also discuss the role of symmetry, which is preserved in our model 
scattering geometry but most likely broken in experimental conditions. 

PACS numbers: 73.63Rt, 73.23.Ad, 73.40.Cg 


I. INTRODUCTION 


Since the observation of zero bias anomalies (ZBA’s) 
for Ce adatoms on silveili^ and Co adatoms on golcP by 
scanning tunneling microscopy (STM), the Kondo effectP 
of magnetic adatoms has become a subject of extremely 
high interest. STM measurements of the Kondo effect 
offer the possibility of achieving exquisite external control 
over a paradigmatic strongly correlated system. Several 
of the established Kondo systems have been reviewed in 
Ref. m Despite the apparent simplicity of these systems, 
a full theoretical description is still lacking. Indeed, great 
efforts using both ab initio and many-body approaches 
have been made, but several open issues still exist. 

The Kondo effect in magnetic adatoms has been suc¬ 
cessfully treated within an Anderson model approaclP 
(mostly with a single impurity or bita l of dz^ symmetry; 
only recently has the whole d shelP^^ been taken into ac¬ 
count), but the role of the tip is still a subject of debate. 
STM measurements have conventionally been performed 
in the tunneling regime where the tip does not affect the 
results. Recent work^SHIIl have looked beyond tunneling 
measurements to explore the contact regime, where the 
geometric details of the tip and its position above the 
adsorbate can affect the Kondo ZB A. Ab initio calcula¬ 
tions including tip-induced perturbations are needed to 
describe, for instance, the observed progression of the 
Kondo temperature as a function of tip height above a 
Co adatom on Cu( 100 jsHin!. 

STM conduc tance is usually calculated in the Tersoff- 
Hamann modeP^^, in which the tip is described by a 
single s orbital and current flows thanks to the coupling 
between the tip and nearby metallic states. However, 


the situation is more complicated in the presence of an 
adsorbate. The interference between tunneling directly 
into the surface and tunneling via the adsorbate causes 
a ZBA in the STM conductance 


G{V) = Gback + AG 


q^ + 2qV/{kBTK)-l 
[g2 + l][{V/kBTK)^ + l] 


( 1 ) 


with a characteristic Fano lineshap^i^, the fingerprint of 
the Kondo effect; here we introduce a l/(g^ -I-1) factor in 
such a way that AG/Gback represents the signal to back¬ 
ground ratio, which is experimentally found to be on the 
order of 10-30%. The parameter q describes the shape of 
the ZBA (<7 = 0 corresponds to a minimum, q = ±oo to 
a maximum, while intermediate values give rise to asym¬ 
metric lineshapes), while Tk, the Kondo temperature, is 
proportional to its width (/cb is the Boltzmann constant). 

The full STM-adatom geometry can alternatively be 
viewed as a nanocontact and the electrical conductance 
calculated within the Landauer-Bfittiker formalisrrfi^. 
But when Kondo correlations are present, the ballistic 
scattering matrix cannot be obtained from density func¬ 
tional theory (DFT) in the standard semi-local approxi¬ 
mations because the latter do not include the many-body 
correlations responsible for ZBA’s. DFT calculations are 
nevertheless indispensable in singling out the relevant ad¬ 
sorption geometries and electronic degrees of freedom. 

To model Kondo ZBA’s from first principles, recent 
workJi^Hi^ have developed a scheme to quantitatively 
join DFT and many-body calculations via an intermedi¬ 
ate Anderson impurity model (AIM). The model param¬ 
eters are determined by matching the mean-field scat¬ 
tering matrix of the AIM to the ballistic scattering ma¬ 
trix of a spin-polarized DFT calculation. Observables 
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are then obtained by solving the AIM with the many- 
body numerical renormalization group (NRG) method. 
This DFT+NRG scheme has successfully predicted the 
Kondo ZBA of nitric oxide adsorbed on the Au(lll) 
surface, albeit underestimating the experimental Kondo 
temperatur^i^. Predictions have also been made for Go 
impurities and vacancies in carbon nanotubeJiSHlSI p)if_ 
ferent approaches that incorporate many-body correla¬ 
tions into first principles calculations thro ugh dynamical 
mean-field theory have also been proposecpl^. 

Two main approaches have been adopted to calculate 
the Fano parameter q in Eq. Q. 

The first is the “two-path model”ISl^ 2 H 26 ]^ which the 
tip is coupled to the adatom d ^2 orbital via hopping 
and to the surface via tpk, where k denotes conduction 
states of the surface and p the states of the tip, giving 
the expression 

_ t-od + y~!i, ^ t2-\-tiVd^G{Q) 

Ek^pk^4d3Gk ^ ^ ^ 

where Vkd are matrix elements between the d orbital and 
the surface and Gk is the Green’s function of the clean 
surface. By considering tpk, tpd and Vk^ to be energy and 
momentum-independent, and calling them respectively 
ti, t 2 and Vd, the right-hand side of Eq. ® is obtained, 
where only the surface Green’s function G(0) at the Eermi 
energy appears. 

The second approach consists in neglect ing t he cou¬ 
pling of the tip to d states, i.e. t 2 = tpd = 0 ,^ 2201 ! gg that 
just the density of states of the metal is probed by STM 
via ti = tpk, assumed for simplicity to be momentum- 
independent; this leads to the Eano parameter 

_ Sk^pk^kdJiGk 3?G(0) 

^ ” Ek^pk^kd^Gk ^ 9G(0)’ ^ ^ 

which is zero for a particle-hole symmetric band if the 
the momentum dependence of the matrix elements can be 
neglected. This second approach is justified for Go/Au by 
the fact that the experimentally observed Eano resonance 
is independent of tip height .1^ If t 2 were not negligible, q 
would be expected to have non-negligible z dependence 
because ti and t 2 generally have different dependence on 
the tip position R = (i?, ip, z) (z being the height above 
the adatom and Ji the lateral displacement in the angular 
direction p). The i? dependence of q has been observed, 
but it is mainly a consequence of probing variations of 
the surface Green’s function at different positions. 

In this paper, we use our DET-I-NRG scheme to study 
how the Kondo temperature and Fano line shape are af¬ 
fected by the location of the tip in two experimentally 
well-characterized cases: single Go impurities on Gu(lOO) 
and Gu(lll) surface^^^^ii^. Experimentally, it is found 
that on the Gu(lOO) surface Go adatoms show a ZBA 
with Tff = 88K and q = 1.1^ upon moving the tip later¬ 
ally q decreases down to O.fP. When the tip approaches 
the Go adatom, Tx increases to 700K and q to ^ 70 in 
one experimenlP^ and T/f=150K and q ^ 2 in anotheJ^. 


On the Gu(lll) surface, Tx = 54K and q = 0.18; upon 
moving the tip laterally, q decreases to ^ 0,*^when the tip 
approaches the adatom, both Tk and q are unaffectecpi^. 

Our main results can be summarized as follows. Spin 
polarized DET calculations show that on both surfaces 
the spin state of Go is S' = 1, each of two magnetic d 
orbitals contributing approximately one Bohr magneton. 
On the (100) surface these orbitals are inequivalent; the 
dz 2 orbital is found to have a much higher Kondo tem¬ 
perature than the d^ 2 _y 2 orbital. The effect of the tip 
is to increase the hybridization of orbital d 2,2 by pushing 
the adatom into the surface, as well as itself providing 
another source of hybridization. On the (111) surface, 
the two magnetic orbitals are degenerate, and their hy¬ 
bridization does not increase as the STM tip gets closer, 
due to symmetry and structural reasons. Gonsequently, 
the Kondo temperature does not vary appreciably all the 
way from the tunneling regime to the contact regime. On 
both surfaces a precise determination of the Kondo tem¬ 
perature, which depends exponentially on the parameters 
of the AIM, is beyond the capabilities of our method. 
Both the Kondo temperature and the Eano parameter q 
are affected by numerous fine details, comprising the elec¬ 
tronic structure of the surface and adatom, the details of 
the tip-adatom-surface nanocontact, in particular, how 
strongly symmetry is broken by the tip, and possibly 
non-equilibrium and multi-orbital (beyond 2) effects. 

The paper is structured as follows. In Section|n]we de¬ 
scribe our method; in Sections |III| and [TV] we present our 
results for Go/Gu(100) and Go/Gu(lll), respectively; 
and in Section 0 we discuss the conclusions of our work. 


II. OVERVIEW OF THE METHOD 

We have employed, with a few simplifications, the 
method presented in Refs. to which we refer for 

further details. 

First we perform a self-consistent, fully relaxed calcu¬ 
lation of the electronic properties of the scattering region 
(as shown in Figs. and by density functional theory. 
This is constituted by a 3 x 3 Gu supercell in the xy plane 
with a Go coverage of 1/9; in the z direction we use 8 Gu 
layers plus a “pyramid” of 5 Gu atoms to simulate the 
STM tip for the (100) surface; for the (111) surface, we 
use 7 layers and a 4-atom pyramid. At this coverage 
the interaction between periodic replicas of the adatom 
is small. The calculations are carried out with the stan¬ 
dard plane-wave package QUANTUM ESPRESSd^ us¬ 
ing the generalized gradient approximation (GGA) to the 
exchange-correlation functional in the parametrization of 
Perdew, Burke and Ernzerhof^. The plane wave cut¬ 
offs are 30 Ry and 300 Ry for the wave functions and 
charge density, respectively. Integration over the two- 
dimensional Brillouin zone is accomplished using a 6 x 6 
grid of k-points and a smearing parameter of 10 mRy. 

After obtaining the self-consistent electronic structure, 
the conductance in the z direction is calculated using 










3 


the PWCOND routine.^ Scattering eigenchannels and 
eigenvalues depend on the transverse momentum 
hence, we restrict our conductance calculations to the 
single most representative \ixy points: B = ^( 5 , 5 ) for 
the (100) surface, and K = J(|,0) for the (111) surface. 
This procedure introduces small systematic errors in the 
estimation of parameters but has the advantage of keep¬ 
ing the computational effort low. We verified that the 
error with respect to a more accurate calculation with 
5x5 \^xy points is less than a few percent. In the above 
expressions, L =7.77A is the length of the supercell in 
the X and y directions, set to three times the equilibrium 
nearest-neighbor distance for bulk Cu for our pseudopo¬ 
tential, 2 . 59 A, slightly larger than the experimental value 
2.56A. 

With the knowledge of the scattering eigenvalues 
it is possible to compute the energy-dependent Fano 
factoil^ 

p _ Sri ^w(f ~ ^n) 

Sn 

which is ex perim entally accessible through noise 
measurements EMII 

In the final step, an AIM is built in such a way that 
it reproduces the DFT scattering phase shifts as closely 
as possible when solved in the Hartree-Fock (HF) ap¬ 
proximation. The Kondo temperature can be estimated 
after the AI M is solved by numerical renormalization 
group (NRG)P^I^. In contrast to Refs. ITHI and IH?! where 
a phase-shift analysis was performed to determine the 
Fano parameter q, here, due to the additional compli¬ 
cation of the dependence on transverse momentum and 
lack of even/odd symmetry along the z direction, the 
lineshape is inferred directly at the DFT level by looking 
at the energy-dependent transmission eigenvalues, and 
fitting them with a Fano lineshape. 

This procedure is repeated for a series of tip-surface 
distances to show how the Kondo temperature and line- 
shape vary in going from the tunneling to contact regime. 

III. CO/CU(100) 

In this section we present results for the Co/Cu(100) 
system, which is found experimentally to have Tx = 88 
K and q = 1.13 in the tunneling regim^. These val¬ 
ues are modified in the contact regime, where Tx grows 
up to 700 K and q up to ^70,^^. Ref. [S] found that q 
does not go beyond 2 - a discrepancy which is probably 
due to the different nanocontact geometries in different 
experimental conditions. 

A. DFT results 

We find that the Co adatom adsorbs in the hollow po¬ 
sition, with 4 nearest-neighbor Cu atoms. In this config¬ 
uration the symmetry group is and Co 3d orbitals 


are split into 3 singlets {dz 2 with symmetry Al, dx 2 -y 2 
with symmetry Bl, dxy with symmetry B2) and a dou¬ 
blet [dxz and dyz with symmetry E); the 4s orbital has 
Al symmetry. In our scattering geometry the tip is built 
so as not to lower the C^y symmetry of the adatom plus 
surface system. We discuss possible consequences of this 
approximation later, since symmetry is not preserved in 
real experiments. 

CCA calculations, both with and without the tip, show 
the presence of two magnetic orbitals, dz 2 and dx 2 -y 2 , in 
agreement with Ref. 1371 while the 4s orbital is highly 
hybridized, almost spin-unpolarized and singly occupied; 
the electronic configuration is thus 3(i®4s^, and the mag¬ 
netic moment is close to 2/1^, in agreement with, for 
example. Ref. |3H1 Table |T] reports some structural and 
electronic data for this system at different tip-surface 
separations. Structural data are in good agreement with 
Ref. [Sni Table |T] covers the approximately-known range 
of experimental tip-surface separations. The tip-surface 
separation could not be further reduced because already 
at the smallest value reported in Table B dtip—sur — 
4.12A, £^2 and Uz 2 are affected by large errors. The 
source of error is the breakdown of our procedure for 
estimating parameters when the occupation of the dj ,2 
orbital grows significantly greater than 1 and part of the 
polarization is transferred to the dxz and dyz orbitals. 
This is expected since with increased coordination Co 
approaches the configuration it has as a bulk impurity, 
where the polarization is shared by all d orbitalpSI. 

Figure shows representative density of states and 
transmission eigenvalues for dup-sur = 5.58A. The inter¬ 
ference between the s and dz 2 orbitals in the down-spin 
Al channel gives a Fano lineshape with g ~ 1 centered at 
an energy around ^ ~ 0.5 eV; all other channels give 
a negligible contribution to the total conductance. 


B. Anderson model 

In building an effective AIM, only the dz 2 and dx 2 -y 2 
magnetic orbitals are retained, each coupled to a linear 
combination of conduction states with the same symme¬ 
try, Al and Bl, respectively. Within this approximation, 
the adatom-surface Hamiltonian is 

ddad—sur ^ ^ T 

k(T 

+ ^ ^ Iki(''kidr'^io- + d^^C\^i„) -\- 

kcr 

-f ^ {e,n, + U,nln\) + 

i=Al,Bl 

-\-JJ\2^z2Bjx 2— y2 JSz2 ' Sx2— y2 ^ (5) 

where we have introduced the single-particle energies Cki 
and fermionic operators Ckio- and associated with 
conduction states with momentum k, symmetry z, and 
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dtip — sur 

1 

0 

0 

dCo—sur 

dco—nn 

m 

€z2 

Uz2 

r.2 

tj.2 y2 

Uz.2 y2 

to 

1 

to 

Ul2 

J 

Tk,z2 

—y 

2 qz2 

9 

F 

7.60 

6.01 

1.59 

2.43 

2.08 

-4.83 

3.06 

0.183 

-4.89 

3.24 

0.147 

1.56 

1.28 

340 

50 

1.20 

0.01 

0.99 

5.58 

3.91 

1.67 

2.47 

2.10 

-4.67 

2.98 

0.165 

-4.80 

3.26 

0.130 

1.50 

1.29 

290 

10 

1.19 

0.31 

0.68 

5.17 

3.41 

1.76 

2.51 

2.10 

-4.54 

2.91 

0.170 

-4.67 

3.23 

0.116 

1.44 

1.30 

410 

3 

0.76 

0.61 

0.53 

4.73 

2.84 

1.89 

2.56 

2.08 

-4.47 

2.86 

0.195 

-4.60 

3.23 

0.101 

1.40 

1.32 

600 

0.5 

0.09 

1.06 

0.42 

4.48 

2.61 

1.87 

2.54 

2.05 

-4.58 

2.88 

0.211 

-4.65 

3.22 

0.104 

1.40 

1.32 

1000 

1 

0.03 

1.27 

0.38 

4.12 

2.42 

1.70 

2.48 

1.96 

-5.47* 

3.36* 

0.232 

-5.23 

3.22 

0.122 

1.63 

1.33 

1100 

25 

0.01 

1.36 

0.34 


TABLE I. Parameters for Co/Cu(100) at different tip-surface separations dup-sur- the distance dup-Co between the tip and 
the Co adatom, dco-sur between the Co adatom and the surface, dco-nn between the Co adatom and its 4 nearest neighbors, 
the total magnetization m of the unit cell in fiB, the on-site energies e ^2 and e^ 2 _y 2 , the Hubbard repulsions U ^2 and U^ 2 _y 2 , 
the hybridizations T ^2 and T^ 2 _y 2 , the inter-orbital Hubbard repulsion U 12 , the Hund exchange J, the calculated Kondo 
temperatures and Tif^^ 2 _y 2 in K, the Fano parameter q^ 2 , the DFT conductance g — G/Go =9^ + 9'^ and Fano factor 

F at the Fermi energy; distances are in A and energies in eV. ‘Here, e ^2 and U ^2 have large errors (see text). 


L = 7.77k 



FIG. 1. Unit cell of the scattering region for Co/Cu(100) for 
dtip — sur ~ 7 . 6 OA (a) in the yz plane and (b) in the xy plane, 
(c) Brillouin zone for the supercell in the xy plane. 


spin O’, the on-site energies Ci and fermionic operators 
dia and associated with impurity orbitals {dAi = 
dsi = dj^ 2 _y 2 ), the hopping elements Vki between a con¬ 
duction state and d orbital, the on-site Hubbard repulsion 
Ui for dz^ and d^-i-y-i orbitals, the inter-orbital Hubbard 
repulsion 17i2 between dz 2 and d^^-y^, and the Hund 
exchange J > 0 between dz 2 and dz: 2 _y 2 . We denote 
< = n, = 1 Hd\^(Ty^d,^, where 

(T is the vector of Pauli matrices cr = CTj,). 


A generalized Tersoff-Hamann model is used for the 
interaction of the tip with the surface and the d orbitals 


on the adatom: 

Htip = + 

per 

+ 'y ^ ^2z(Cpjo.dicr + dj^Cpio-) + 
per 

+ ^ ki^Cpia), (6) 

pker 

where Cpi and Cpio- denote the single-particle energies and 
destruction operators associated with conduction states 
of the tip with momentum p, symmetry i, and spin a. 
The usual approach only takes into account the A1 sym¬ 
metry channel, i.e. the tip-surface hopping tiAi and t 2 Al^ 
which is motivated by the fact the apex atom of the tip 
(Cu in this case) has a single s orbital at the Fermi energy. 
In contrast, we allow for conductance through channels 
with different symmetry. Since our geometry preserves 
the symmetry of the Co d states, the total conductance 
gtot = GtotjGq^ where Gq = jh is the quantum of con¬ 
ductance, is the sum of the conductance from A1 and HI 
channels (we ignore all other symmetry channels as they 
do not contribute to the ZBA): 

gtot = 9i — gAi + 931- (J) 

i=Al,A2,Bl,B2,E 

We are aware that symmetry is not preserved in real ex¬ 
periments, because the tip cannot be expected to have 
the ideal pyramid-like shape we have assumed; how¬ 
ever, it is a good starting point for studying the prob¬ 
lem. Should symmetry be broken, channels with different 
symmetries would start interfering, leading to a modified 
lineshape. However, we believe the disagreement we find 
with the experimentally determined q is only be partially 
due to this approximation: since a single Kondo temper¬ 
ature is relevant (see next subsection), the main effect of 
interference would be to slightly modify the hopping pa¬ 
rameters from the tip to the surface and the adsorbate, 
only weakly affecting our estimate of q. 
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FIG. 2. Co/Cu(100): (a) Spin-polarized density of states 

at Co 3d and 4s atomic orbitals constructed from scattering 
states at ka;,j, = B = |;(|, |) (positive DOS=spin up, nega¬ 
tive DOS=spin down); (b) and components of the DPT 
transmission, transmission eigenvalues and 

ti = t'*') and the shot noise Fano factor. Energies are with 
respect to the Fermi energy. 


We now introduce the hybridization functions rf(e) 
due to the coupling with the surface, r-(e) due to the 


coupling with 

the tip. 

and the total ri(e): 



-E 

k 

^(e - eki)Ek^i ^ 

(8) 

r‘(e) 

-E 


(9) 

r.(e) 


"brj(e), i = dz2, dx2—y2, 

(10) 

where the expressions 

TrpsVi^ and T^Ptt^i are valid 

as long 


as we assume energy independent quantities {psi and pu 
denote the density of states of symmetry i at the Fermi 
energy for the surface and the tip, respectively). Our 
method only allows us to infer the total ri(e), which we 
assume to be energy-independent, and simply call F^. 
Also, we assume that the density of states of conduction 


electrons is flat, and extends from —Di to Dp. 




g(A-|£|) 

2A 


i = A1,B1, 


( 11 ) 


and we take Dai = Dbi = 7eV as estimated from the 
density of states at the Fermi energy of s electrons for 
Cu atoms on the clean surface. 

In principle, the 4s orbital of the Co atom should also 
be taken into account: 

'y ( “b y ( l^Als(’^CT^kAlfT “b ^kAlcr'^rr) “b 

kcr 

~t~ts y ] (Cpylicr®cr + sl-CpAlcr) ~ 
pAlo" 

-Ss■Y.Js^S,, (12) 

i 

with on-site energy e^, destruction operator So- for spin tr, 
hybridization matrix elements V^ais with A1 conduction 
states, coupling tg to A1 states of the tip, and exchange 
coupling Jsi with di states. However, due do its large on 
(Fs ^ 2eV), the 4s orbital can be taken as part of the 
A1 conduction band, thus effectively enhancing tiAi in 
Eq.i for the dz 2 orbital. 

In practice, since we are mainly interested in the Kondo 
temperature, instead of solving Eqs. and together 
by NRG, we always solve Eq. alone but replace the 
hybridization F® due to the surface with the total hy¬ 
bridization Fi = F® -b F* due to the surface plus tip. As 
explained in Sec. |TTj the lineshape is approximated dur¬ 
ing the DFT step without resorting to the model Hamil¬ 
tonian in Eqs. and 

The parameters in Eq. [^are then fixed by trying to re¬ 
produce as closely as possible the GGA results within the 
HF approximation of the AIM in the wide-band limitfi^, 
which gives: 

e] = ei + Uiuj + ^ U.jUj - ^ Jijrrij/d, (13) 

j 3 

ej = ei + Uiu] + ^ ^ JijTrij/d, (14) 

3 3 

where we sum over all j ^ i atomic orbitals (nf, with 
a =t, is the occupation of orbital i in the spin chan¬ 
nel cr; Tii = n\ -\- n\] rrii = n} — nj). The linewidths 
Fi are taken from the down-spin density of states by fit¬ 
ting each impurity orbital with a Lorentzian; this actu¬ 
ally gives ri(e^), but we assume ri(e^) ~ ri(0). J is 
assumed to be constant in the d-shell. It is inferred from 
the energy splitting of the dxy orbital (which has a very 
low magnetization mxy O.Od/is), induced by the total 
magnetization m via exchange interactions: 


J = 


2(A - ft 

\^xy 


m 


(15) 


Hubbard repulsions Ui are taken from the splitting of 
magnetic orbital, once J is known: 

ej - ej = Uimi +^{m-TTii), i = d^i.dx^-y^- (16) 
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The inter-orbital Hubbard repulsion U 12 is approximated 

bjOSl 


C /12 =C/ave- (17) 

where Have is the average of the Hubbard repulsion over 
the ^^2 and dx 2 -y 2 orbitals. Finally, the on-site energies 
ei are fixed from the orbital energies ef together with the 
knowledge of Ui and C/ 12 : 

ej + et Ui • ■ j j ■ _L ■ 

62 — ^ ; ‘^1 J — ^2^ , dx‘^ — y^ •) J ^ 

(18) 

Here, the numerical values of ef, nf, and rrii are taken 
from DFT. The results of this procedure is reported in 
Table H 


C. Kondo temperature 

In order to estimate the Kondo temperature, Eq. [^for 
the dj ,2 and orbitals is solved by a two-channel 

NRG calculation. Each orbital, with on-site energy e^ 
and Hubbard repulsion Ui, is coupled to its own Wil¬ 
son chain through the full broadening F^ = Ff -|- F*, 
which takes into account the interaction with both the 
surface and tip; the two channels are coupled by Hund 
exchange coupling J and inter-orbital Hubbard repulsion 
C/ 12 . Tke Kondo temperature is obtained by computing 
the spectral function for both impurity levels, and taking 
the halfwidth of the zero energy resonance. 

GGA predicts the down-spin orbitals to lie just above 
the Fermi energy (see Fig. [^. When translated into an 
AIM, this means that ej ^ Ci + Ui ^ 0, which leads to 
high particle-hole asymmetry. As a consequence, NRG 
predicts the magnetic orbitals to be almost in the mixed- 
valence region, which explains why the values of Tx in 
Table |T] are higher than the experimental ones. However, 
we believe this high particle-hole asymmetry is a spurious 
effect, which could be amended by resorting to some more 
sophisticated method, such as GGA-|-ljS]. For example, 
we found that a small value of C/ ^ 0.5 eV in the GGA-I-U 
approach is enough to reproduce the experimental Tk in 
the tunneling regime. 

Our value of Fj ,2 ~ 0.18 eV in the tunneling regime is 
comparable, but somewhat lower than other values found 
in the literature: 0.20 eV in Ref. [28land 0.24 eV in Ref. El 
Even though our method of computing the hybridization 
from the broadening of scattering states is in principle 
more accurate than a simple estimation from the den¬ 
sity of states after a self consistent calculation, since it 
involves the interaction of impurity levels with a contin¬ 
uum of states, hence requiring no artificial broadening, 
the approximation of using a single k^, ^ point can easily 
lead to a ~ 10% error, as a consequence of the interaction 
among periodic replicas of the impurity in the x - y plane, 
which could be alternatively reduced by the use of larger 
supercells. In addition, the way the nuclear relaxation is 


performed—either taking into account magnetic effects 
or ignoring them—is found to be another source of un¬ 
certainty. For example, performing a spin-unpolarized re¬ 
laxation increases from 0.183 to 0.208 eV with respect 
to the standard spin-polarized calculations we use in this 
paper. This differs with earlier worlP^, where magnetic 
and nonmagnetic calculations were found to yield similar 
relaxed atomic coordinates due to a cancellation between 
adatom-substrate and adatom-tip interactions; see also 
Ref. [TTl For comparison, we find the relaxed Go-surface 
distance without the tip to be 1.59A when magnetism is 
taken into account, and 1.49A otherwise; in Ref. EH this 
value was found to be 1.5lA, regardless of magnetism. 

As such, a precise evaluation of the Kondo temperature 
remains a challenge. However, we stress that we obtain 
the correct growth of in passing from the tunneling 
regime to the contact regime, as reported in Table |T] and 
Fig.0 In addition, in the intermediate regime, i.e. for 
g ^ 0.3, Tx^z^ decreases slightly, which is probably true 
experimentally (see Fig. 3a of Ref. [Tn|l . This is a con¬ 
sequence of the fact that the Co adatom is attracted by 
the tip, and therefore pulled away from the surface (see 
Table m column dco-sur', this agrees well with Ref. EH). 
This causes a decrease in the hybridization r 22 , which 
is compensated only in the contact regime, when the tip 
is close enough to the Co atom to be considered as an 
additional nearest-neighbor atom. At this point, rj ,2 and 
Tx^z 2 are greatly enhanced. However, this is only true for 
the fij .2 orbital, which is probably the one producing the 
experimental ZBA. In the case of dx 2 _y 2 , the tip greatly 
decreases Tx^x^-y^ when approaching the adatom, be¬ 
cause the decreased hybridization of the orbital with the 
surface due to the increased surface-adatom distance is 
not compensated by an additional hybridization with the 
tip, due to symmetry mismatch; should symmetry be 
broken, things would not be different, since the dx 2 -y 2 
orbital lies flat on the surface, thus hardly coupling with 
the STM tip no matter how this approaches the adatom. 
This is true until the tip really “pushes” the atom closer 
to the surface, eventually enhancing Tx^x^-y^ too. We 
stress that, in contrast to Refs. El and uni we attribute 
the change of the Kondo temperature mainly to changes 
in the hybridization, rather than in the on-site energy 
and Hubbard repulsion. 


D. Lineshape 


As remarked, in our symmetry-preserving geometry, 
the total lineshape is the sum of the A1 and B1 con¬ 
ductances in Eq. However, after solving Eq. with 
the parameters shown in Table III we find that Tx^z^ ^ 
Tx,x 2 -y 2 ', moreover, gAi ^ gsi oecause the dx^-y^ or¬ 
bital is flat and only couples to the second layer of the 
tip for symmetry reasons. This suggests that most of the 
experimental Kondo signal is due to the (i ^2 orbital and, 
in what follows, we shall assume gtot = gAi and ignore 
(/Bi, together with all other symmetry channels, which 
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do not carry ZBA’s and contribute very little to the con¬ 
ductance. In any case, we find that the Fano parameter 
associated with the orbital dx 2 _y 2 should always be much 
larger than 1, and this would show up as an additional 
weak anomaly superimposed to the standard one. 

Applying the phase-shift analysis described in Ref. Uni 
turned out to be too cumbersome for this system, due to 
its intrinsic 3D character, so instead we estimated q from 
the shape of the energy-dependent DFT transmission 
coefficient, which shows interference at energies around 

^ 0.5 eV. If we assume that the hopping parameters 
in Eq. are weakly energy dependent, the DFT line- 
shape can be a good approximation to the ZBA, which 
involves interference at the Fermi energy. Unfortunately, 
when we do so, the agreement with experiment is not 
always good. In the tunneling regme, we obtain g ^ 1, 
which nicely matches experiment^. However, when go¬ 
ing into the contact regime, we find a decrease of q. In 
this regime the conductance G is close to the unitary 
value Go, thus the interference between the s and or¬ 
bitals can only be destructive, leading to a dip in the con¬ 
ductance, while in experiments the opposite is observed 
(q increases slightlji^or strongljff^. 

A possible reason for this disagreement in the contact 
regime is that the junction is formed in a different way 
than we have modeled it. Other sources of error in our 
calculations are the inclusion of non-equilibrium effects 
and our assumption of energy-independent parameters; 
see also the conclusions in Sec. |Vl 


IV. CO/CU(lll) 

In this section we report our results for Co on the 
Cu(lll) surface, which has Tk = 54 K and q = 0.18 
in the tunneling regim^. These values remain almost 
unchanged in passing to the contact regim^^^l. 


A. DFT results 

Once again, we find that Co adsorbs in the hollow po¬ 
sition, this time with 3 nearest-neighbor Cu atoms in a 
configuration with symmetry. The d orbitals split 
into a singlet with symmetry Ai {dz 2 ) and two doublets 
with symmetry E {dai^ 2 j <^^ 1 , 2 ); the s orbital has Al 
symmetry. The electronic configuration is 3d®4s^, and 
the total magnetic moment is close to 2/rs, just like on 
the (100) surface. We again model the tip to preserve 
the symmetry (C^y). 

If we take Cartesian axes as in Fig. the doublets can 


be written in the following way: 

dal = cosddxz +sm9dxy, (19) 

da2 = cosOdyz + s\ii9dx2-y2, ( 20 ) 

dpi = — sin 9dxz + cos 9dxy, (21) 

dp 2 = — sin 9dyz + cos 9dx2-y2 , (22) 


where the dai,pi orbitals are odd with respect to the sym¬ 
metry operator Px : x ^ —x, while the da 2 ,p 2 orbitals 
are even. From DFT calculations, it turns out that the 
dal ,2 doublet is magnetic, while (i/ 31,2 is fully occupied; 
moreover we find 0 = 0.70 rad. 

In Table m we show some structural and electronic 
data for different tip-surface separations. In Fig. |^we 
show density of states and transmission eigenvalues for 
the shortest distance, dup-sur = 4.33A. At K, the de¬ 
generacy between dai and da 2 is weakly broken, so two 
different peaks appear in the density of states and in 
transmission eigenvalues at ~ —0.5eV {dp orbitals) and 
~ 0.5eV {da orbitals) for down-spin electrons in the E 
channels. 



FIG. 3. Unit cell of the scattering region for Co/Cu(lll) for 
dtip — sur — 7. 81A (a) in the yz plane and (b) in the xy plane, 
(c) Brillouin zone for the snpercell in the xy plane. 


B. Anderson model 

The full atomic Hamiltonian is 

H = ^ {tim + Ui rii-^ nil) - 

i—z^ ,al,a2,pi,P2 

^ ( Jij T ^ ( UijniTlj j (^ii) 

j<i j<i 

with on-site energies tai = £02 = e/ 3 i = e /32 = e/ 3 , 
Hubbard repulsion Uai = Ua 2 = Ua, and Upi = Up 2 = 
Up, inter-orbital Hubbard repulsion Uij and Hund ex¬ 
change coupling Jij. After dropping fully-occupied or¬ 
bitals, which means keeping only dai and da 2 , we intro- 
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dtip — sur 

dtip —Co 

dco — sur 

dco—nn 

m 


Ua 

Va 

Gi2 

J 

Tk 

Q 

g 

F 

7.81 

6.01 

1.80 

2.40 

2.20 

-5.20 

3.26 

0.164 

1.90 

1.09 

70 

7.1 

0.01 

0.99 

5.81 

3.92 

1.89 

2.44 

2.20 

-5.33 

3.34 

0.168 

1.98 

1.09 

50 

6.4 

0.47 

0.53 

4.33 

2.44 

1.89 

2.44 

2.19 

-5.39 

3.28 

0.157 

1.92 

1.09 

100 

4.5 

1.06 

0.10 


TABLE II. Parameters for Co/Cu(lll) at different tip-surface separations dup-sur'- the distance dup-co between the tip and 
the Co adatom, dco-sur between the Co adatom and the surface, dco-nn between the Co adatom and its 3 nearest neighbors, 
the total magnetization m of the unit cell in units of ps, the on-site energy ea, the Hubbard repulsion Ua, the hybridization 
Pq,, the inter-orbital Hubbard repulsion U 12 , the Hund exchange J, the calculated Kondo temperature Tk in K, the Fano 
parameter q, the DFT conductance g = and the Fano factor F at the Fermi energy; distances are in A and energies in 

eV. 


duce metallic states, to get 


Had 


— sur 


ekC^j^Ckicr + 

kfj 

+ y ' + 

kcr 

i—(y.l,Oi2 

+ yZ {(^ani + Uan\n\) + 

2=0:1,0:2 

H“t^l2^al^o2 J *^ol ' ^0.2- 


In addition to the d orbitals, the s orbital is half-hlled 
and highly hybridized, exactly as on the (100) surface, 
with the same Hamiltonian (Eq. 12), but it is irrelevant 
when dealing with Kondo physics: it only contributes to 
the conductance in the Al channel which shows no ZBA, 
and does not interfere with the da orbitals, which have 
different symmetry E. 

Once again, a generalized Tersoff-Hamann model 


Htip — ^ + (25) 

per 

2=o1,o2 

y ' + '"kicr'^pi'r) + (26) 

pker 

2=o1,o2 

^ ^ ^2i{Cpi(j^io' "t“ d\^Cpicr)f (27) 

per 

2=o1,o2 

is used and the hybridization functions 

r«(e) = 7r^^(e- ek)14a ^ ^PsV^ ( 28 ) 

k 

rL(e) = 7ry^i5(e-ek)tL ^ (29) 

k 

Taie) = Kie) + riie) (30) 


are introduced and approximated as energy-independent 
quantities F^, r(,, r^. 

The parameters are fixed in the same way as for the 
Co/Cu(100) case; we get Da = 5 eV from the density of 
states of surface s orbitals; J is hxed from the splitting of 
P orbitals, which have mp = 0.06. Numerical values are 
reported in Table [III The value Tq, ^ 0.16 eV is slightly 
below the value 0.18 eV of Ref. H31 


C. Kondo temperature 


The Kondo temperature is obtained by solving Eq. 
with NRG for the model parameters in Table [ll| We used 
the total hybridization Tq, including contributions from 
both the surface and tip. 

Since the dai and da 2 orbitals are degenerate, there is 
a single Kondo temperature as reported in Table [lT| Spin- 
orbit effects will lift this degeneracy, leading to two differ¬ 
ent Kondo temperatures, but the effect should be small. 
As for the previous case, we overestimate the Kondo tem¬ 
perature due to the excessive particle-hole asymmetry 
which comes from GGA, but in this case the disagree¬ 
ment is not too bad. 

When the tip is brought down to the surface, we 
find that the Kondo temperature first decreases slightly 
(even though the difference is below the accuracy of our 
method), because the adatom is pulled farther away from 
the surface, then weakly increases. The s orbital of the 
tip is of different symmetry than the a orbitals, so T^ = 0 
to a first approximation. The effect of the tip-induced re¬ 
laxation of the Co adatom and its neighbors is found to 
be negligible, in contrast to the (100) case, so T* is ba¬ 
sically unaffected by the tip. This is in good agreement 
with a similar analysis in Ref. [11] and with experimen¬ 
tal results that show a constant Tx as a function of the 
tip positioiJSl. At a distance dup-sur = 4.33A, the Co 
atom is pushed towards the surface, and at this point the 
Kondo temperature starts to increase, but this regime is 
probably not reached in experiments. 

In Fig. we show the evolution of the conductance g 
and Kondo temperature Tk^ as a function of the tip- 
surface distance for both Co/Cu(100) and Co/Cu(lll). 
The reported values of are to be seen mostly as up¬ 
per estimates, since the use of GGA-I-U, as mentioned, 
would decrease the Kondo temperature. In fact, our val¬ 
ues overestimate the experimental Kondo temperature. 
It must be stressed that most of the uncertainties are 
systematic, so they affect all the data in the same direc¬ 
tion. 
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FIG. 4. Co/Cu(lll): (a) spin-polarized density of states 
at Co 3d and 4s atomic orbitals constructed from scattering 
states at = ^(|,0) (positive DOS means spin up, negative 
DOS spin down); (b) and t'*' components of the DFT trans¬ 
mission, transmission eigenvalues tJ and t)" (eigenvalues close 
to 1 correspond to the A1 channel while two lower eigenval¬ 
ues are from the E channels), and the shot noise Fano factor. 
Energies are with respect to the Fermi energy. 


D. Lineshape 

When the tip is above the adatom, symmetry is pre¬ 
served, and the total conductance is 

gtot = <?i — 2 ^ 0 , + Gai- (31) 

i—Al,a,j3 

In the Al channel, there is no ZBA, because only the s 
orbital is involved, the ( 1^2 orbital being completely filled. 
In the E channel, there is a ZBA due to da orbitals, but 
the signal should be small because they do not couple to 
the s orbital of the tip. DFT predicts that these channels 
give a peak in the conductance (g I, see Fig. [^, the 
coupling of the tip to the orbitals being much higher than 
to Cu states, according to the two-path model, but still 
much lower than the coupling to the Co s orbital. How¬ 


ever, when compared to experiment, where a dip {q ^ 0) 
in conductance is observed, with a strong signal, this is 
wrong. Moreover, on the basis of the DFT results, one 
would expect to see considerable changes in the lineshape 
when moving the tip laterally in the xy plane, due to 
symmetry breaking, but this is not seen either. Finally, 
when the tip approaches the surface, the lineshape is un¬ 
affected, remaining a minimum, again in contradiction 
with what one would expect from the DFT results. 

This probably means that symmetry is unimportant, 
most likely because the tip breaks it. This implies that 
one can observe E orbitals even when the tip is above the 
adatom, and when moving the tip laterally there is no 
symmetry breaking and no significant change in the sig¬ 
nal. However, this is not sufficient to explain the exper¬ 
imental results because one would still expect the ZBA 
to be a peak. Another interpretation might be that since 
the magnetic orbitals have E symmetry, while tunnel¬ 
ing only happens through metal states of Al symmetry, 
these acquire a dip in their density of states, leading to 
q ^ 0, as suggested in Ref. ISH however, metal states 
of Al symmetry should carry no, or at most very weak, 
Hondo signal. 

Instead, it is likely that surface states are responsible 
for the dip. According to DFT, the tip “sees” the mag¬ 
netic orbital and not the metallic states, while it should 
be the other way around. This is simply because the Co 
adatom is closer to the tip than to the Cu surface atoms 
(on the (100) surface things are different, because the 
Co s orbital, which is effectively part of the conduction 
band, can interfere with the magnetic orbital.) However, 
surface states, if taken into account, might prove to be 
more prominent than the d orbitals, leading to a dip in 
conductance. This interpretation is also suggested by the 
“quantum mirage”!^ experiment, which shows how sur¬ 
face states can carry Hondo information even far from 
the impurity, with the same lineshape as when the tip 
is above the adatom. Our DFT slab calculation can in 
principle include surface states, but to describe them cor¬ 
rectly we would need a much larger supercell. 

The role of surface states has been discussed in the 
literature. The total hybridization from the surface 
r® = Tsurf + ^buik is the sum of the hybridization 
from surfa ce states and bulk states. While it is gener¬ 
ally agreecPl2lIlIIl2HlII^ the only exception being Ref. H51 
based on calculations and experimental hints (for exam¬ 
ple the lack of appreciable changes in the Hondo temper¬ 
ature at step edges and defect^, where surface states are 
deeply affected), that Fbu/fc is much larger than Fgu^y, by 
up to a factor of one hundrecP^, not much is known about 
the relative magnitude of tip-surface and tip-d orbitals 
coupling, which is what controls the lineshape. It is only 
known that for a clean surface about two thirds of the 
current flows into surface stated. Ref. 021 argues that 
surface states can give an important contribution to the 
conductance even in the presence of adsorbates. 

In any case, it must be stressed that the usual assump¬ 
tion that the magnetic orbital is of dj ;2 character is found 
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to be wrong in this case. This makes the usual Tersoff- 
Hamann approach fail, because due to symmetry mis¬ 
match there can be no conductance from the s-state of 
the tip into the magnetic orbital, or to the linear combi¬ 
nation of conduction states to which the magnetic orbital 
is coupled, which leads to the paradox that when the tip 
is directly above the adatom, it should give almost no 
signal. 

V. CONCLUSIONS 

We have demonstrated the application of our 
DFT-I-NRG methocP^ to Co adatoms on Cu surfaces 
probed by an STM tip, trying to describe the lineshape 
and the Kondo temperature both in the tunneling and 
contact regimes. 

The calculations show the severe difficulty in predict¬ 
ing the details of experimental Kondo anomalies, espe¬ 
cially as far as the lineshape is concerned, while the trend 
in the Kondo temperature, if not its absolute value, is 
reliable. Several different issues might cause these dis¬ 
crepancies. 

First of all, our description of the tip is surely oversim¬ 
plified. Yet, experimental measurements are only weakly 
dependent on the choice of the tip, at least in the tun¬ 
neling regime, so in principle this should not be a large 
source of error. 

We believe that on the Cu(lOO) surface the essential 
remaining point is to correct the estimation of the param¬ 
eters which control the lineshape: the asymmetric ZBA 
in the tunneling regime can be understood in terms of 
the interference between the s and ^ 2,2 orbitals of the Co 
adatom. Issues may arise in the contact regime due to 
the geometrical details of the contact, which may differ 
from our model. The momentum dependence of the hop¬ 
ping parameters such as 14,, Cipk, ^ 2 ip, which we have 
ignored, might also play a role, as well as non-equilibrium 
effects. Also, in the contact regime, which is not far from 
a bulk impurity situation, all d orbitals start to become 
magnetized, thus actively entering conduction processes, 
and making our two-orbital model insufficient. 

On the Cu(lll) surface, in contrast, more work should 


be done to include surface states. This might reverse the 
sign of the ZBA (giving a dip instead of a peak). Also, 
one should seek to understand if the E symmetry of the 
magnetic orbital, instead of A1 as usually assumed, can 
affect the lineshape, as suggested e.g. in Ref. |37l 

On both surfaces our estimate of the Kondo temper¬ 
ature would be improved by correcting the excessive 
particle-hole asymmetry brought about by plain GCA, 
for example by using GGA-fU^. Moreover, our many- 
body model may be improved. For example, one could 
take into account spin-orbit effects, correlated hopping, 
double hopping, and other two-body interactions, or keep 
all d orbitals, considering that some of them are almost 
but not completely filled, especially in the contact regime. 
In addition, the energy dependence of the hybridization 
functions ri(e) is likely to have some impact. Finally, we 
note that the co-existence of many magnetic solutions 
(not only of lowest energy presented in this paper) real¬ 
ized at different Co adatoms (on both Cu surfaces ) could 
be responsible for statistical spread in Kondo tempera¬ 
tures and lineshapes. 

We emphasize that, according to our GCA results, Co 
has spin S' = I on both Cu surfaces, showing that the 
usual assumption of S = 1/2 is probably wrong. How¬ 
ever, this mistake might not have a big impact on the final 
result if one of the two magnetic orbital has a much lower 
Kondo temperature than the other one, as on the (100) 
surface, or if the two magnetic orbitals are degenerate, as 
on the (III) surface. Of course, a detailed quantitative 
approach cannot overlook this fact. 

In any case, despite several works claiming that the 
Kondo physics of adatoms is fully understood, we believe 
that further effort is needed to completely understand, 
or at least describe satisfactorily from first principles, 
the Kondo behavior of Co adatoms on Cu surfaces, and, 
more generally, of magnetic adatoms on metallic surfaces. 
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